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ABSTRACT 

We present new UVES spectra of a sample of 15 cool unevolved stars with and without 
detected planetary companions. Together with previous determinations, we study Be depiction 
and possible differences in Be abundances between both groups of stars. We obtain a final sample 
of 89 and 40 stars with and without planets, respectively, which covers a wide range of effective 
temperatures, from 4700 K to 6400 K, and includes several cool dwarf stars for the first time. 

We determine Be abundances for these stars and find that for most of them (the coolest 
ones) the Bell resonance lines are often undetectable, implying significant Be depletion. While 
for hot stars Be abundances are aproximately constant, with a slight fall as T e g decreases and 
the Li-Be gap around 6300 K, we find a steep drop of Be content as T e g decreases for T c ff < 
5500 K, confirming the results of previous papers. Therefore, for these stars there is an unknown 
mechanism destroying Be that is not reflected in current models of Be depletion. 

Moreover, this strong Be depletion in cool objects takes place for all the stars regardless of 
the presence of planets, thus, the effect of extra Li depletion in solar-type stars with planets when 
compared with stars without detected planets does not seem to be present for Be, although the 
number of stars at those temperatures is still small to reach a final conclusion. 



Subject headings: stars: abundances - stars: fundamental parameters 
satellites: formation - stars: atmospheres 



planetary systems - planets and 



1. Introduction 

Light elements are important tracers of stel- 
lar internal mixing and rotation. Since they are 
burned at relatively low temperatures they con- 
strain how the material inside stars is mixed with 
the hotter interior. Rotation and angular mo- 
mentum loss are among the leading processes to 
explain the mixing that leads to depletion of light 
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elements in solar-ty pe stars (e.g. IStephens et al 
19971 : lBouvierll2008l) although gravitational waves 
may also affect the abundances of those elements 
( Montalban fc Schatzman .2000). However, avail- 
able models for evolution of Be (considering ro- 
tation) do not predict a significant depletion 
of Be during the main sequ ence for stars with 
6000 K > T Gff > 4000 K (jPinsonneault et al. 
1990h . On the other hand, models which take 
into account gravitational waves predict signifi- 
ca nt Be depletion only for stars co oler than 4500 
K (|Montalban fc Schatzmanll2000l) . 



In a recent work, llsraelian et al. (l2009h 



firmed that Li was severely depleted in solar-type 
stars (with T e ff between 5600 K and 5850 K) with 
planets when compared with similar stars without 
detecte d planets although th is result is contro- 



versial ( Baumann et al. 201fj Sousa et al. 2010l) 



for a complete discussion see lDelgado Mena et al 
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(|201ll ). This difference in Li abundance seems 
to be related to the different rotational history 
of both groups of stars d ue to the pre sence of 
protoplanetary disks (e.g. Bouvier 20081 ). How- 
ever, beryllium needs a greater temperature to be 
burned so we would expect to see the onset of this 
effect in cooler stars where convective envelopes 
are deep enough to reach those higher tempera- 
tures. 



In a previous paper. iDelgado Mena et al.l(|201lh 
found two cool planet host stars with an extra de- 
pletion of Be when compared with analog stars 
without detected planets. This encourages us to 
try to investigate this process in cool stars. In 
this work we continue that analysis by extending 
the sample with 15 new cool stars. We refer the 
reader to that paper for further information and 
a more extensive introduction. 



2. Observations and spectral synthesis 

In this study we obtained high resolution spec- 
tra for 15 new stars with magnitudes V between 6 
and 10 using the UVES spectrograph at the 8.2- 
m Kueyen VLT (UT2) telescope (run ID 86.D- 
0082A) between October 2010 and March 2011. 
The dichroic mirror was used to obtain also red 
spectra and the slit width was 0.5 arcsec. These 
new spectra have a spectral resolution R ~ 70000 
and S/N ratios between 100 and 200. All the data 
were reduced with the pipeline of UVES/VLT. 
Standard background correction, flat-field and ex- 
traction procedures were used. The wavelength 
calibration was made using a ThAr lamp spec- 
trum taken during the same night. Finally we 
manually normalized the continuum by dividing 
the spectra by a spline function with three pieces 
(and the parameter low — reject set to 1) in the 
whole blue region (3040A-3800A). This normaliza- 
tion does not present any bias for stars with and 
without planets. When plotting together observed 
spectra of stars with and without planets we only 
had to multiply the flux of comparison stars by 
0.9-1.1 in order to make them match up. We note 
that in our previous works we have analyzed high 
metallicity (up to 0.4) solar type stars with Teff 
down to 5300 K for which we could make a good 
normalization. The main problem is the normal- 



ization of spectra of stars with Teff less than ~ 
5300 K. These stars could have spots and inho- 
mogeneous atmospheres that make the spectral 
synthesis more difficult. 

The uniform stell ar atmospheri c parameters 
were taken from ISousa et al.l (2008) with typical 
errors of 25 K for T e g, 0.04 dex for logg, 0.03 km 
s _1 for £ t and 0.02 dex for metallicity. We refer to 
that work for further details in these parameters 
and their uncertainties. 

Be abundances were derived by fitting the spec- 
tral region around the Be II line at 3131.06 A 
and then using the Be II line 3130.42 A to check 
the consistency of the fit. We used an empiri- 
cal line list from Garcia Lopez et all (|l995 ) tuned 
to reproduce the solar spectrum (see Table []}. 
These synthetic spectra were convolved with a ro- 
tational profile. We made a standard LTE analysis 
with the revised ve rsion of the sp ectral synthesis 
code MOOG2002 (jSnedenl Ewm and a grid of 
Kurucz ATLAS9 at mospheres with overshooting 
(jKurucz et al.lll993h . Examples of synthetic spec- 



tra and the parameters used in the synthesis are 
shown in Figure [TJ 

The final sample is composed of 70 and 30 
stars with and without plan et s , respectively, from 



Santos et all (|2002l l2004al0 ); iGalvez-Ortiz et al 



(2011), 14 stars with planets from lDelgado Mena et al 
(|2011l ) and 5 and 10 stars with and without de- 
tected planets, respectively, from this work. This 
gives a total sample of 89 stars with planets and 
40 comparison sample stars. All Be abundances 
for these 89+40 stars were analyzed by our team 
using the same methodology, making this a very 
uniform sample. 
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Fig. 1. — Synthetic spectra (red lines) and observed spectra (dots) for the planet host star HD 93083 and the 
stars without detected planets HD 213042, HD 35854, HD 15337, HD 8389A and HD 21019. The position 
of Be lines are indicated by the arrows. 
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Fig. 2. — Panel a: Spectral synthesis of HD 213042 using different line lists. Panel b: Spectral synthesis of 
HD 16270 showing the contributions of atomic, molecular and Be lines. Panel c: Spectral synthesis of HD 
63454 with solar line list and different values of Be abundance. Panel d: Spectral synthesis of the star Rl in 
the young cluster IC 2602 using different line lists. 
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Table 1 — Continued 
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0.592E-01 


3129.652 


41.1 


1.321 


0.115E+00 


3132 


288 


25.0 


4.332 


0.316E+00 


3129.763 


40.1 


0.039 


0.331E+00 


3132 


355 


23.0 


1.043 


0.984E-04 


3129.774 


24.0 


2.708 


0.344E-02 


3132 


392 


108.0 


1.990 


0.131E-03 


3129.857 


24.0 


2.968 


0.845E-02 


3132 


405 


25.0 


3.373 


0.818E-02 
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3. Analysis 

In general, Be abundances for stars cooler than 
5200 K are probably not reliable since in this 
regime Be lines are barely sensitive to changes in 
the abundance. In Figure [1] we can observe that 
for the coolest stars the fits are not good. At those 
temperatures Mn I line at 3129.037 A dominates 
the feature and the presence of Be is negligible (see 



also IGarcia Lopez et al.l 119951 ) . We make several 



tests to try to improve those fits. 

In Figure [5] we show several spectral synthe- 
ses. In panel a we have used different line lists 
for HD 213042. The orange dashed line is a fit 
made with original Kurucz line lisld. It is clear 
that the Mn-Be feature cannot be well reproduce 
with the original value of log gf for Mn line even 
if we increase Be abundance (light blue dashed- 
pointed line). Furthermore, the lambda of the 
whole feature do not match Be line, so this star 
do not present so much Be and anot her line is re- 



quired to fit the observed spectrum. iPrimas et al 



(|19971 ) used an artificial Fe line at 3131.043 which 
in our case does not help to fit the spectrum (pur- 
ple dashed-three pointed line) since this line gets 
stronger in metallic stars. Another option is to 
increase the gf of Mn line at 3131.037 A This 
was first proposed by IGarcia Lopez et "ah ( 1995 ) 
to reproduce solar spectrum and we also used this 
modified line in our previous works on Be. The 
green line represents a fit with Kurucz line list 
and this modified line. We note that for our syn- 
theses we have use d Mn measured abundances by 
Neves et al.l (|2009f ) with the same models. This fit 
is better but the feature is still stronger than ob- 
served in its red wing although the Be abundance 
used in the syntehsis is negligible. The synthetic 
spectra of these cool stars present strong molecu- 
lar lines. If we decrease O and C abundances (red 
line) we can get a better fit though this does not 
affect Be line at 3131.06 A. 



core shifted towards higher A which would not fit 
the observed spectrum. We can again get a better 
fit if we decrease O and C abundances (red line). 
However, as the stars get cooler we need to use 
lower C and O abundances, ~ -0.6 dex for O and 
~ -0.9 dex for C in this star. This is not a realis- 
tic approach but fortunately C and O abundances 
hardly affect the Be feature. 

In panel c we show three synthe ses for HD 

63454 with the solar line list derived bv lGarcia Lopez et al 
(|l995l ) and decreasing C and O abundances -0.3 
and -0.2 dex respectively. The fits with log Be = 
0.6 and 0.3 dex presents a stronger line than ob- 
served, so we can adopt an upper limit of 0.0 dex 
for this star. Certainly, Be line is not very sensitive 
at these low temperatures and we think that we 
may be overestimating Be depletion. However, it 
is imposible to fit the spectra using high Be abun- 
dances so we can put an upper limit in Be content 
though we cannot calculate accurate abundances. 
Moreover, the sensitivity of Be line does not seem 
to be related to T e ff since some cool young objets 
prese nt strong Be lines with abundances similar t o 
solar (jSmilianic et al.ll201l l ; iRandich et al.ll2007l) . 
This is the case of the star Rl in the young cluster 
IC 2602 of 46 M yr and solar metallicit y. This star 
was observed bv lSmiljanic et al.l ( 20111 ) who found 
log Be = 1.25. We have analyzed this star using 
the same spectrum and thr ee different li n e lists 
(see panel d of Figure [5]): IPrimas et al. ( 1997 ) 
line list, used by these authors, gives log Be = 
1.22; our solar line list gives log Be = 1.37 and 
Kurucz line list with Mn line modified and de- 
creasing C and O abundances by 0.3 dex gives log 
Be = 1.32, all of them in perfect agreement with 
the value previously found. Therefore, this test 
probably suggests that our line list does work for 
cool stars which present a measurable quantity of 
Be, and our fittings should be valid at least to put 
an upper limit in Be abundances. 



The different contributions of molecular an 
atomic lines is shown in panel b for star HD 16270. 
The Mn-Be feature is practically filled with atomic 
lines, that is, Mn line. Therefore we can put an 
upper limit in Be abundance since increasing Be 
abundance would result in a stronger line with a 



2 http:/ /kurucz. harvard.edu/line lists.html 



4. Discussion 

In this section we will discuss the different 
trends of Be abundances with the effective temper- 
ature, the metallicity and the oxygen abundance. 
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Table 1 — Continued 



A(A) 


Atomic number 


X (eV) 


gf 


A(A) 




Atomic number 


X (eV) 


9f 


o 1 on no a 

3129.934 


OA 1 

39.1 


O A 1 A 

3.414 


0.955E+01 


01 no 

3132 


517 


68.1 


1 A AO 

1.402 


C\ OOAT7 1 i A1 

0.320E+01 


3129.937 


108.0 


1.609 


1.955E-02 


3132 


518 


26.0 


7.168 


0.316E+00 


3129.943 


73.0 


0.697 


0.724E-01 


3132 


532 


24.1 


6.805 


0.387E-03 


3129.968 


64.1 


1.172 


0.627E+00 


3132 


579 


26.1 


7.495 


0.161E-02 


3129.974 


90.1 


1.287 


0.172E+00 


3132 


583 


108.0 


1.612 


0.314E-01 


3130.056 


40.0 


0.519 


0.200E+00 


3132 


591 


58.1 


0.295 


0.288E+00 


3130.063 


26.1 


13.018 


0.416E-02 


3132 


594 


42.0 


0.000 


0.237E+01 


3130.075 


108.0 


2.295 


0.0029174 


3132 


596 


23.1 


2.900 


0.859E-01 


3130.126 


108.0 


0.842 


0.794E-02 


3132 


631 


607.0 


0.510 


0.210E-04 


3130.145 


108.0 


1.987 


0.0141579 


3132 


656 


73.0 


0.491 


0.110E+00 


3130.157 


22.0 


5.941 


0.346E+00 


3132 


657 


27.0 


2.878 


0.240E-04 


3130.202 


25.1 


4.801 


0.193E+00 


3132 


710 


22.0 


5.954 


0.171E+00 


3130.254 


106.0 


0.521 


0.0660693 


3132 


725 


25.1 


6.177 


0.215E-02 


3130.257 


23.1 


0.348 


0.513E+00 


3132 


788 


25.0 


4.268 


0.316E+00 


3130.281 


108.0 


0.250 


0.134E-01 


3132 


794 


8.2 


36.895 


0.933E+00 


3130.290 


106.0 


0.521 


0.0731139 


3132 


809 


23.1 


2.510 


0.297E-01 


3130.340 


58.1 


0.529 


0.705E+00 


3132 


816 


108.0 


1.947 


0.313E-03 


3130.353 


27.1 


2.985 


0.465E-03 


3132 


822 


24.0 


3.122 


0.322E+00 


3130.370 


106.0 


0.033 


0.111E-01 


3132 


845 


108.0 


1.947 


0.348E-02 


3130.376 


22.0 


1.430 


0.275E-01 


3132 


864 


28.1 


2.865 


0.223E-03 


3130.407 


108.0 


1.756 


0.324E-03 


3132 


865 


108.0 


0.686 


0.578E-03 


3130.408 


108.0 


1.670 


0.103E-03 


3132 


878 


44.0 


1.317 


0.174E+00 



*ine list used bv iGarcia Lopez et al. ( 1995I ) and in our previous works. The dissociation energy for OH, CH, 
NH and CN molecules is 4.39, 3.46, 3.47 and 7.65 eV respectively. 
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Table 2 
Stars analyzed in this work. 



Star 


T e ff 


log g 


6 


[Fe/H] 


V 


planet 


Spectral type a 


log e(Be) 


log e(Li) 




[K] 


[cm s~ 2 ] 


[km s- 1 ] 














HD2638 


5198 


4.43 


0.74 


0.12 


9.44 


yes 


G5 


0.49 


<0.16 


HD8326 


4971 


4.48 


0.81 


0.02 


8.70 


no 


K2V 


i— 0.16 


<0.09 


HD8389A 


5283 


4.37 


1.06 


0.34 


7.84 


no 


KOVCN+2 


0.16 


<0.73 


HD9796 


5179 


4.38 


0.66 


-0.25 


8.81 


no 


K0V 


0.27 


<0.17 


HD11964A b 


5332 


3.90 


0.99 


0.08 


6.42 


yes 


G9VCN+1 


0.55 


1.41 


HD15337 


5179 


4.39 


0.70 


0.06 


9.10 


no 


K1V 


0.58 


<0.42 


HD16270 


4786 


4.39 


0.84 


0.06 


8.37 


no 


K3.5Vk: 


<-0.32 


<0.03 


HD21019 b 


5468 


3.93 


1.05 


-0.45 


6.20 


no 


G2V 


0.22 


1.39 


HD27894 


4952 


4.39 


0.78 


0.20 


9.42 


yes 


K2V 


<-0.38 


<0.22 


HD35854 


4928 


4.46 


0.54 


-0.13 


7.74 


no 


K2V 


<-0.31 


<-0.22 


HD40105 b 


5137 


3.85 


0.97 


0.06 


6.52 


no 


K1IV-V 


<-0.12 


1.40 


HD44573 


5071 


4.48 


0.80 


-0.07 


8.46 


no 


K2.5Vk: 


0.75 


<-0.01 


HD63454 


4840 


4.30 


0.81 


0.06 


9.37 


yes 


K3Vk: 


<-0.32 


<-0.03 


HD93083 


5105 


4.43 


0.94 


0.09 


8.33 


yes 


K2IV-V 


<-0.14 


<0.16 


HD213042 


4831 


4.38 


0.82 


0.08 


7.66 


no 


K5V 


<-0.40 


<0.06 



a Values taken from Simbad 
b Evolved stars 
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Fig. 3. — Observed spectra and difference in fluxes for six pairs of planet-host stars (red lines) and stars 
without detected planets (blue dashed lines). The position of Be lines are indicated by the arrows. 



4.1. Is Be depleted in stars with planets? 

We have seen in the previous section that ab- 
solute Be abundances determined with spectral 
synthesis are probably not reliable for cool stars, 
preventing a comparison between stars with and 
without detected planets. In order to search for 
differences in Be abundances between those stars 



we proposed to c ompare directly their spect ra in 
a previous paper (jDelgado Mena et al.ll201lh . We 
selected pairs of stars with similar stellar param- 
eters so if there were a difference between their 
spectra in the Be region it should be due to a 
difference in Be abundance. In that paper we 
presented two planet-host stars, HD 330075 and 
HD 13445, with an extra Be depiction when cran- 



io 



v a ^ 





< 


5150 K • 


5150 < T e , 


< 


5600 K • 


5600 < T ef 


< 


5900 K • 


5900 < T„ 


< 


6500 K • 



Log N(Li) 



Fig. 4. — Left panel: Be abundances as a function of effective temperature for dwarf stars with and without 
detected planets from this work (red filled and blue open triangles, respectivel y) and dwarf stars with (re d 



filled circles) and without planets (blue open circl es) from previous studies (jSantos et all I2002 , 2004alT< 
Galvez-Ortiz et al. 2011; Dclga do Mena et al. I l201ll) . The three evolved stars from this work are depicted 



wi th squares. The Sun is de noted by the usual symbol. The dashed lines represent 4 Be depletion models 
of Pinsonneault et al. ( 199fjf) (Case A) with different initial angular momentum for s olar metallicit y and a n 
age of 1.7 Gyr. The solid line represents an assumed initial Be abundance of 1.26 Santos et al.l ( 2004c ). 
The dotted line repre sents the Be depletion iso c hrone for 4.6 Gyr taken from the models including mixing 
by internal waves of Montalban fc Schatzman ( 2000f) . Right panel: Be abundances as a function of Li 
abundances. Filled and open circles are stars with and without planets, respectively, from previous surveys. 
Filled and open triangles are stars with and without planets, respectively, from this work. Filled and open 
squares are evolved stars with and without detected planets, respectively, with measured Be abundance in 
this work. Colors denote different temperature ranges. 



pared to analogous stars without detected planets, 
showing that the effect observed in Li abundances 
for solar-type stars with giant planets might also 
occur for Be in cooler stars. 

In this work we present new spectra for 15 cool 
stars (see Table[2]) with and without detected plan- 
ets. U sing these stars and ot hers from previous 
works dSantos et all l2004allo!: iGalvez-Ortiz et al. 
bOllUDelgado Mena et al.ll201n) ~ we made 13 new 
pairs of analogous stars with differences in T e ff, 
log g and [Fe/H] lower than 50K, 0.15 dex and 0.10 
dex respectively. We also note that v sini values 
of these stars are very similar. Some examples of 
these couples are shown in Figure [3] We can see 
in these plots that stars with and without planets 
present similar Be features and the differences in 
flux between spectra are very low around Be fea- 
tures. In the previous paper we found three stars 
from the comparison sample with clearly higher Be 
abundance than the planet-host star HD 330075, 



whose abundance is log e(Be) < -0.25. Now we 
have two more stars with very similar parameters. 
HD 44573 has a higher Be abundance (see Table 
[2]) and the difference in fluxes between both stars 
is around 5<r in the position of the two Be lines. 
However, HD 8326 has a spectrum very similar to 
HD 330075 and Be abundances are almost equal. 
Therefore, planet host stars and comparison stars 
can have similar Be abundances. We note, how- 
ever, that our sample of cool stars is small and 
although unlikely, it might be possible that we are 
only observing the fraction of comparison stars 
with a strong Be depletion. A similar situation 
is seen for solar-type stars without detected plan- 
ets, where 50% have high Li abundances while the 
other 50% do have its Li depleted like planet-hosts. 

4.2. Be versus T G ff 

In left panel of Figure |4] we plot the derived Be 
abundances as a function of effective temperature 
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• T eff < 


5 


50 


K 






• 5150 


< 


T eir 


< 


5600 


K 


• 5600 


< 


T efr 


< 


5900 


K 


• 5900 


< 


T ef( 


< 


6500 


K 



[Fe/H] 



5600 K < T„„ < 



log Be = 1.14 ■ 
51 50 K log Be = 1.13 • 
log Be = 1.25 



0.07 

0.38 



F H" 
Fe/H = 
Fe/H = 



[Fe/H] 



Fig. 5. — Left panel: Be abundances as a function of metallicity. Filled and open circles are dwarf stars with 
and without planets, respectively, from previous surveys. Filled and open triangles are dwarf stars with and 
without planets, respectively, from this work. Colors denote different temperature ranges. Right panel: Be 
abundances as a function of metallicity for dwarf stars with 5600 K < T e ff < 6150 K and detections of Be. 
We have overplotted two linear least square fittings for planet-host stars (red filled circles, red line) and stars 
without detected planets (blue open circles, blue dashed line). 



for planet-host stars in our sample (see Table[5]) to- 
gether with previous samples. In this plot we have 
removed subgiants and giants (except the three 
analyzed in this work) to avoid evolutionary effects 
in the abundances. In this selection we took the 



spectr al types of the stars from lGalvez-Ortiz et al 
(|201lh . With this new sample of stars we have 



completed the coolest part of the plot and now we 
can see the behaviour of Be abundances in a wide 
range of effective temperatures. 

As mentioned in previous papers, the Be abun- 
dances decrease from a maximum near T Q g = 
6100 K towards higher and lower temperatures, 
in a similar way as Li abundances behave. In the 
high temperature domain, the steep decrease with 
increasing temperatures resembles the well know n 
Be gap for F stars (e.g. iBoesgaard fe Kindl2002h . 
The decrease of the Be content towards lower tem- 
peratures is smoother and may show evidence for 
continuous Be burning during the main sequence 
evolution of these stars. 

When we move to the coolest stars, we observe 
a strong Be depletion regardless of the presence 
of planets. Hence, the planetary formation pro- 
cesses that affect Li depletion in solar anal og stars 
( Israelian et al. 20091 ; Gonzalez et alj|2010l) do not 
seem to take effect in the Be abundances. 



On the other hand, we have found several 
stars with undetectable or very weak Be lines. 
The steep decrease of Be abundances for T c ff 
< 5500 K is in contradiction w i th m odels of 
Be depletion (jPinsonneault et al. 1990), which 



predict either constant or increasing Be abun- 
dances as T c ff decreases (see Figure H}. Even 
taking into account mix i ng by internal waves 
( Montalban fc Schatzmanl I2000T ) . Be depletion 
is still lower than predicted. Be abundances 
have also been studied in clusters of different 
ages. 



Randich et al.1 (|2007l ) found that stars in 



the 150 Myr old cluster NGC 2516 with 5000 
< T c ff < 5500 K have not depleted Be while 
Hyades stars (600 Myr) of the same T e ff present 
some Be depletion, something that cannot be 
expla ined by models wit h only convective mix- 



ing. iGarcia Lopez et al.l (|1995I ) also found evi 



dence of Be depletion for two cool Hyades mem- 
bers. Although uncertainties in Be abundandes 
for the coolest stars are large and it might be 
some kind of systematic effect due to the Mn I 
A 3131.037 A line, many of them seem to have 
their Be totally destroyed. Therefore, none of 
those models including convective mixing, rotation 
or gravitational waves seems to fit the observed 
Be abundances, at least at these cool tempera- 
tures. Other depletion mechanisms have been pro- 
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posed to e xplain Li destruct i on, such as tachoc line 
diffusion (|Brun et al.l Il999t IPiau et all 120031 ) or 
episodic a ccretion of planetary mat erial during 
the PMS (|Baraffe fe Chabrierl [2O10h though Be 
has not been analyzed. However, accretion of 



tational hi stories dSantos et al.l 2004c ). In a re- 



cent work, iTakeda et al.l ([2011!) found four stars 



planetary material in the ZA MS (jTheado et al 



2010l : lTheado fc Vauclairlboilh can lead to an en- 
hanced depletion of Li and Be, as well. These 
models can explain the difference in the Li abun- 
dances between stars with and without planets 
and might also explain the observed Be differ- 
ences in a pair of cool planet-hosts when com- 
pared to several stars wit hout detected planets 
( Delgado Mena et al.ll201ll ). On the other hand 
the depletion produced by episodic accretion does 
not seem to be enough to explain the strong Be 
depletion observed in cool stars and the presence 
of planetary material would be necessary in all 
stars though many of them do not have detected 
planets. 



4.3. Be versus Li 

A beryllium versus lithium diagram can give 
us information about the depletion rates in main- 
sequence stars. In right panel of Figure [4] this 
relation is shown for all the stars in our sample. 
In general Be abundances increase with increasing 
Li abundances and stars with and without plan- 
ets (filled and open symbols, respectively) behave 
in a similar way. If we take into account the ef- 
fective temperatures we can divide the sample in 
four groups. The first one, with the hottest stars 
(light blue points), present both high Be and Li 
abundances except for two stars with T e ff around 
6300K which fall in the Be-Li dip region, where 
the depletion of those elements is atributted to 
slow mixing and depends on age and temperatur e 
(jBoesgaard fc KingM2002l : lBoesgaard et al.ll2004al ) . 



though one of them, HD 120136 presents a nor- 
mal Li abundance (log e(Li) = 2.04). In the sec- 
ond group (orange points), with effective temper- 
atures between 5600 and 5900 K, Be abundances 
remain high although Li abundances present dif- 
ferent depletion rates. There are several objects 
with low Be abundances and solar temperatures 
(specially the comparison star HD 20766 T c ff = 
5733 K, log e(Be) < -0.09), which could form 
some kind of "Be-gap" where the depletion of Be 
may be related to different pre-main sequence ro- 



with solar temperatures and Be line undetectable. 
Those stars are also depleted in Li and have low 
v sini values indicating that slow rotation can be 
the cause for that strong Be depletion. 

In the third group (blue points), composed by 
stars with 5150K < T G ff < 5600K, Be abundances 
begin to decrease and Li is severely depleted ex- 
cept for two objects, HD 11964A and HD 21019, 
with particularly high Li abundances considering 
its temperature. These two objects are evolved 
stars, therefore a possible explanation for the high 
Li content is that they have just left the main 
sequence, where they were hotter than they are 
now, but there has not been time for their Li 
and Be to become strongly depleted. Other pos- 
sible explanations are a dredge up effect from a 
"buffer" below the former main sequen ce convec- 



tive envelope ( Delivannis et al. 1990h 



or accre- 
tion of planetary metal-rich material, although 
recent models indicate that accretion of plane- 
tary material might de stroy Li instead of p r oduc- 
ing an enhanceme nt (jBaraffe fc Chabrierl 1201ft 
Theado et aTll2010l) . Finally, the coolest objects 
of our sample (red symbols) are depleted in both 
Be and Li, although some stars have preserved 
some Be. There are two objects with anomalous 
high Li content. HD 74576 seems to be a young 
star, something that proba bly justifies its high 
Li content (see discussion in ISantos et al] l2004d l 
while HD 40105 is an evolved star that might 
have suffered some of the processes proposed for 
HD 11964A and HD 21019. We note that we have 
removed from the plots all the evolved stars from 
previous works. Only the three subgiants observed 
for this work are analyzed here but t here are sim- 



ilar cases in the whole sample (see ISantos et al 
l2004d : iGalvez-Ortiz et aT1l201ll ). 



This fig ure confirms what found inlSantos et al 



(|2004al lch; IGalvez-Ortiz et al.l (|201ll) . Be and Li 
burning seems to follow the same trend. Both 
elements increase their content as the tempera- 
ture rises, although Li depletion begins at 5900 
K while strong Be depletion starts below 5500 K. 
Stars with and without planets behave in a similar 
way. 
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4.4. Be versus [Fe/H] 

In the left panel of Figure [5] we show Be abun- 
dances as a function of metallicity for different 
temperature ranges. We can see that the stars 
are equally distributed regardless of their temper- 
ature and for higher metallicities the dispersion 
in Be abundances increases, mainly due to objects 
with 5150 K < T cff < 5600 K. It is well known that 
Be abundances increase with metalli city for [Fe/H] 
< -1 with a steep slope nea r 1 (e.g. iRebolo et al 



fore, oxygen abundances can provide us with com- 
plementary information about the galactic evo- 
lution of Be. We c o mpile oxygen abundances 



1988 : Boesgaard et al.ll2009h . but for higher metal- 
licities this relation is not so well defined. 

Boesgaard et al. (2004b) found a slope of 38 



for stars with m etallicity between -0 .6 and 0.2, in 
agreement with iTakeda et al. ( 2011 ) whose slope 
value is 0.49, although they take into account only 
solar analog s with -0.3 < [Fe/H] < 0.3. On the 
other hand, iBoesgaard et al.l (|2009h argued that 
stars of one solar mass with solar metallicities 
(their most metallic star has [Fe/H] = 0.11) match 
the slope of 0.86 for metal-poor stars. 

In the right panel of Figure [5] all the dwarf stars 
of our sample with 5600 K < T g- < 6150 K are 
plotted. We have chosen this limit in temperature 
in order to remove most of the upper limits in Be 
abundances and those stars from the Li-Be dip. 
We have also removed stars from the solar Be-gap. 
We find a slope of 0.13 for planet-host stars and 
0.07 for comparison sample stars. These values are 
lower than previously found by other authors but 
this could be possibly due to the lack of stars with 
[Fe/H] < -0.2 in our sample in comparison with 
the high number of metal-rich objects. It seems 
that stars with solar and supersolar metallicities 
do not follow the same trend as metal-poor stars, 
and Be has been produced at a lower rate as the 
Galaxy has evolved. However, we note that Be 
abundances present higher uncertainties in metal- 
rich stars that might be affecting these trends. If 
we have a look in the upper envelope, which might 
be more consistent with low er metallicity stars 
(around [Fe/H] = -0.5), we obtain a fit with a slope 
of 0.38, more consistent with previous results. 

4.5. Be versus [O/H] 

Beryllium is produced by spallation reactions 
between galactic cosmic rays (GCRs) and the 
CNO nuclei in the interstellar medium (see e.g. 
Tan et al. 20091 and references therein). There- 



from Ecuvillon et al 



Galvez-Ortiz et al 



d2006l) (see the final lists 



20111 ). derived with the OH 



bands in the near-UV since this region is closed 
to Be lines. We use the same temperature ranges 
than in previous section, 5600 K < T e g < 6150 K. 
Thus, none of the stars analyzed in this work are 
included in this plot because all of them are cooler 
than 5600 K. How ever, 6 planet-host stars f rom 
our previous work (jDelgado Mena et al.ll201ll ) are 
hotter and we have displayed them as green cir- 
cles in Figure [6] Oxygen abundances for these 
stars have bee n obta ined in the same way as 
Ecuvillon et al. (l2006h and are presented in Ta- 



ble n 

In the middle panel of Figure [6j Be abundances 
are plotted as a function of [O/H]. We can see 
that the slopes of planet hosts and stars without 
detected planets are slightly different, possibly due 
to the low number of comparison sample stars in 
this range. We have removed the star HD154857 
since its O abundance is the lowest of this sample 
and its content in Be is considerable lower than 
other stars of similar oxygen abundance, hence it 
would increase the slope by 0.3. We make again 
a fit of the upper envelope for stars with plan- 
ets to take into account the lack of stars with low 
[O/H] in our plot. These trends might be influ- 
enced by the effect of T c ff in Be abundances. To 
remove this effect, we have made a new Be vs. 
[O /H] plot (lower panel of Figure [6]) with "cor- 
rected" Be abundances, calculated as log e(Be)- 
log e(Be(T e ff)), where log e(Be(T ff)) is obtained 
from a linear fit of Be as a function o f tempera- 
ture in the same T c ff range (upper panel of Figure 
[6|). However, the slope is still the same for stars 
with planets, it seems that there remains a corre- 
lation between Be and O abundances, which may 
indicate that the sensitivity of Be to the Teff is 
not affecting at all this correlation. On the other 
hand, the comparison sample is so small at this 
temperature range, that no strong 

conclusion can be extracted from the different 
slopes obtained before and after the Be- Teff linear 
correction. 

Nevertheless, these slopes are considerable 
lower than previously found for metal-poor stars, 
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Table 3 

Oxygen abundances for stars with T cff > 5400 K from [Delgado Mena et a"l~I (|201ll ). 



Star 


[O/H]! 


[0/H] 2 


[0/H] 3 


[0/H] 4 


[0 /ii] final 


log e(Be) 


HD2039 


0.17 


0.12 


0.17 


0.29 


0.18 


1.44 


HD4203 


0.08 


0.01 


0.22 


0.20 


0.13 


1.18 


HD73526 






0.20 


0.25 


0.22 


1.24 


HD76700 


0.16 


0.16 


0.28 


0.35 


0.24 


1.08 


HD154857 


-0.23 


-0.20 


-0.18 


-0.18 


-0.20 


0.54 


HD208847 


0.13 


0.08 


0.05 


0.10 


0.09 


1.28 


HD216770 


-0.01 


0.04 


0.06 


0.06 


0.04 


0.66 



even t he slope of the upper en velope in the middle 
panel. iBoesgaard et al. ( 19991 ) obtained a slope of 
1.45, using; stars with -3 < [Fe/H] < 0.05, while 
Tan et al.l (|2009 ) reported a slope of 1.49, for stars 
with -2.3 < [Fe/H] < -0.6. Therefore, the slope 
for metal-rich stars is much flatter than for metal- 
poor stars as it seems to happen for Be versus 
[Fe/H]. A slope like this would not be in agree- 
ment with models of Be production t hat predict 



quad ratic or linear relations (see e.g. I Tan et al 



2009, and references therein) . We note again that 



oxygen abundances present a high dispersion in 
metal-rich stars and even at low metallicities, the 
use of different O indicators change this trend. 

5. Conclusions 

We present new high-resolution UVES/VLT 
near-UV spectra of 15 stars with and without 
planets in order to find possible differences in Be 
abundances between them and confirm our previ- 
ous results that suggested a greater depletion of 
Be in planet-host stars when compared with stars 
without detected planets. Be needs higher tem- 
peratures than Li to be destroyed so we have to 
search for these differences in cooler stars, which 
have deeper convective envelopes and are able to 
carry the material towards Be-burning layers. We 
have made new couples of analog stars with and 
without planets with the purpose of comparing 
directly their spectra to see possible differences 
in Be abundances. Although in a previous work 
we found a pair of planet-host stars with their 
Be severely depleted, now we have not observed 
important differences in Be abundances between 
both groups of stars. Thus, the effect caused by 



protoplanetary disks and rotational history on ex- 
tra Li depletion in solar-type stars with planets is 
not taking effect in Be abundances, apparently. 

Furthermore, Be abundances in cool stars (T e g 
< 5200 K) are much lower than predicted by the 
models. We have analyzed for the first time a 
considerable number of cool objects and found a 
strong destruction of Be in most of the stars. Al- 
though we cannot provide absolute abundances 
we can give upper values for Be. This gives a 
steep drop of Be abundances as T e ff diminishes in 
contradiction with current models of Be depletion. 

Finally, the slopes of Be abundances as a func- 
tion of [Fe/H] and [O/H] are considerably lower 
than found for metal-poor stars, indicating that 
the production rate of Be may have diminished 
with the evolution of the Galaxy. 
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• • • cC 



stars with planets slope = 0,34 
comparison stars slope = 0.12 
stars with planets upper envelope = 0.52 



* 



stars with planets slope = G. 33 
comparison stars slope = 0.01 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



Fig. 6. — Upper panel: Be abundances as a func- 
tion of temperature for dwarf stars with effective 
temperature between 5600 K and 6150 K with (red 
hlled circles) and without planets (blue open cir- 
cles). Green circles are stars with planets with 
oxygen abundances measured in this work. Middle 
panel: Be abundances as a function of [O/H] us- 
ing O abundances from OH ( Ecuvillon et alj [2006) 
lines for stars with effective temperature between 
5400 K and 6300K. Symbols like in top panel. 
Lower panel: Be abundances corrected of temper- 
ature effect for the same stars than middle panel. 
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